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We address the issue of precisely estimating small parameters encoded in a general linear trans¬ 
formation of the modes of a bosonic quantum field. Such Bogoliubov transformations frequently 
appear in the context of quantum optics. We provide a set of instructions for computing the quan¬ 
tum Fisher information for arbitrary pure initial states. We show that the maximally achievable 
precision of estimation is inversely proportional to the squared average particle number and that 
such Heisenberg scaling requires non-classical, but not necessarily entangled states. Our method 
further allows us to quantify losses in precision arising from being able to monitor only finitely many 
modes, for which we identify a lower bound. 
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I. INTRODUCTION 

Quantum metrology exploits distinctive quantum fea¬ 
tures, such as entanglement, to enhance the estimation 
precision of parameters governing the dynamical evolu¬ 
tion of the probe systems beyond that achievable by clas¬ 
sical means. This enhancement is manifested in the form 
of a scaling gap in precision with respect to the available 
resources (the number of probe systems or the average 
input energy) between the corresponding optimal quan¬ 
tum and classical strategies, and depends on the spe¬ 
cific encoding of the parameter in the Hamiltonian de¬ 
scribing the evolution. In the case where the parameter 
of interest is a multiplicative factor of a local Hamilto¬ 
nian acting on N probes, the optimal quantum strat¬ 
egy provides a quadratic scaling gap in N, known as the 
Heisenberg limit, over the best classical strategy [1, 2]. 
This quadratic improvement is also present for a class 
of quasilocal Hamiltonians [3, 4]. For Hamiltonians in¬ 
volving highly non-local interactions, super-Heisenberg 
scaling is also possible [5-7]. Nonetheless, the paradig¬ 
matic example for a scaling gap is the estimation of a 
phase acquired in one arm of a Mach-Zehnder interfer¬ 
ometer [1], There, the two input modes are subject to a 
particular Gaussian transformation, i.e., a combination 
of beam splitters and a phase shifters. 

Here we are interested in determining the ultimate pre¬ 
cision limits for a more general type of parameter esti¬ 
mation task, namely, where the parameter of interest is 
encoded nontrivially in an arbitrary Gaussian transfor¬ 
mation of not two but possibly infinitely many modes 
of a quantum field. This problem is of broad inter¬ 
est since such transformations, often cast in the form 
of Bogoliubov transformations, feature in a large variety 
of physical systems throughout quantum optics [8] and 


* nicolai.friis@uibk.ac.at 


condensed matter physics [9]. In exchange for allowing 
for a broader class of Gaussian transformations, we re¬ 
strict our approach to a regime of small parameters to 
gain analytical insights. Besides applications such as es¬ 
timating (weak) single-mode or multimode squeezing in 
optical or superconducting (see, e.g., [10, 11]) systems, 
perturbative Bogoliubov transformations of this kind are 
of particular interest for the description of quantum ef¬ 
fects in curved spacetime, such as the Unruh effect [12] 
and the dynamical Casimir effect [13], or analogous re¬ 
alizations thereof [14-18]. A paradigm that highlights 
the challenges encountered in the context of estimating 
relativistic quantum effects [19—23] is the estimation of 
the acceleration of a nonuniformly moving cavity [24]. 
While the parameter in question is small, the Gaussian 
transformation couples all pairs of modes in nontrivial 
ways [25, 26]. Due to the resulting notoriously cumber¬ 
some perturbative calculations, the only known bounds 
on precision involve Gaussian input states. Consequently, 
the ultimate limits on how precisely one can determine 
small accelerations, as well as whether distinctive quan¬ 
tum features provide an improvement, are not known. 

Here we show that, within such a perturbative ap¬ 
proach, there exists a quadratic scaling gap between the 
optimal quantum and classical strategies for a fixed av¬ 
erage input energy. We construct an optimal estimation 
strategy utilizing separable input states of fixed parti¬ 
cle number and boson counting measurements, where we 
use the quantum Fisher information (QFI) as a figure of 
merit for the estimation precision. We provide simple for¬ 
mulas for the leading-order contributions to the QFI for 
arbitrary pure states, both when all modes can be con¬ 
trolled and when only part of the spectrum is accessible, 
and show that, within the considered regime, Heisenberg 
scaling is the ultimate precision limit. Moreover, we iden¬ 
tify the family of states that exhibit Heisenberg scaling 
and show that, while some of these states may be en¬ 
tangled, the crucial feature is their nonclassicality rather 
than their correlations. Furthermore, we lower bound the 
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loss in precision due to tracing out inaccessible modes and 
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provide criteria for minimizing such losses. 
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In its most general form, our metrological protocol can 
be described as follows: We wish to estimate as precisely 
as possible a single real parameter 9 that is encoded in a 
unitary transformation U ( 6 ) acting on the Fock space 
of an arbitrary number of bosonic modes. Our probe 
system is a set of (in principle infinitely many) noninter¬ 
acting harmonic oscillators with creation and annihila¬ 
tion operators a' n and a n , respectively, where n = 1,2,... 
can be a multilabel distinguishing frequencies, polariza¬ 
tions, or other degrees of freedom. The operators sat¬ 
isfy the usual commutation relations [a m ,aj)] = S mn 
and [a m ,a n ] = 0 and the ground state is annihilated by 
all a n , i.e., a n |0) = 0 V n. Arbitrary pure states | if) 
can be decomposed into superpositions of Fock states, 
e.g., | TOfcj ) | nk 2 ) • • • | Pk N }, with fixed numbers of exci¬ 
tations. These Fock states form a basis of the total 
Hilbert space and can be obtained from the vacuum by 
applying the appropriate creation operators, i.e., \nk) = 

(4)7v^|o>. 

For a given input state p and a suitable measurement 
whose outcomes are used to estimate 6 , the precision of 
the estimation of 9, quantified by the variance A 0 of 
the corresponding (unbiased) estimator, is lower bounded 
by the inverse of the QFI X(p(9)) [27] via the (quan¬ 
tum) Cramer-Rao inequality [28, 29] A 9 > l/i Jvl(p(9)), 
where v is the number of repetitions. Intuitively, the QFI 
quantifies how well small changes of the parameter in 
question may be inferred from measurements of the final 
state p{9) after the dynamical evolution of the probes. 
One may express I(p(9)) as 


l{p(9)) 


1 - 

lim 8- 

0 


y/F(p(9),p(9 + d0)) 
d9 2 


(1) 


where F(p(0),p(0+d8)) = (Tr yj\fp(9) p(9 + d9)^p(8)) 2 
is the Uhlmann fidelity between the states p{9) and 
p(9 + dO). If p(9) = \ip(8))(ip(9)\ one simply recovers 
F(p(9),p(9 + d9)) = | (ip(9) | %jj(8 + d9) )| 2 . The latter ex¬ 
pression applies for the case where all modes can be con¬ 
trolled and measured. Expanding \ip(9 + d9)) in powers 
of d9 , one obtains 


FIG. 1. (Color online) Parameter estimation scheme. To 

estimate the parameter 6, a set of modes k = {fci, k 2 , ..., fcjv} 
is prepared in some initial states | ip) k , while all other modes 
are left in the vacuum. The unitary transformation U(6) en¬ 
codes the parameter in the final state | ip ), but measurements 
to determine 9 can only be performed on the marginal pk, 
while information in other modes is lost. 


III. QFI IN THE PERTURBATIVE REGIME 


A. Parameter estimation for pure states 

We now consider cases where 9 is very small and 
close to zero. For this case the transformation U{9) is 
amenable to a perturbative approach and has a Taylor- 
Maclaurin expansion of the form 

U(8) = U w + 9 U w + 9 2 U (2) + 0(8 3 ) , (3) 

where U {0) is a unitary that encodes the free time evolu¬ 
tion of the noninteracting bosonic modes, while for * > 0 
the U (i) represent (nonunitary) perturbations. Conse¬ 
quently, the final state can be written as 

\i>) = |# 0) ) + 9 |^M) + 9 2 \^) + 0(9 3 ), (4) 

where |^ (i) ) = U (i) \ip) (i = 0,1,2,...). Note that |i p) 
and | are normalized, but the vectors | (i > 0) 

generally are not. Substituting (4) into the expression 
for the QFI in (2), we arrive at 

X(|$» = 4(<$< 1 i|$< 1 >> - |(# 0) |^ (1> )| 2 ) +0(8). (5) 

Hence, given that all the bosonic modes are accessible, 
the QFI can be straightforwardly computed for any initial 
(pure) state from just its linear perturbations. 

As all our results utilize the perturbative approach, a 
few important remarks concerning the applicability of 
the perturbative approach are in order. In the latter 
we assume that all higher-order terms in 9 are small 
and are therefore neglected. To ensure such reasoning 
is justified, we regard our approach as valid so long 
as the perturbation to the state | 'i/> (0) ) remains small, 
such that 1 — ^(l ^ (0) ), | ip)) <C 1. This implies that 
9 2 I(\ip))/4 <C 1 and all optimizations will be performed 
with this constraint in mind. 


where | tp(8)) = \if(9)). When some of the modes are 

not accessible and are traced out, Eq. (2) provides an up¬ 
per bound on the QFI. To allow for more specific state¬ 
ments about the QFI, additional information about the 
initial states or transformations encoding 9 is required. 


B. Tracing losses 

Thus far we have assumed that all bosonic modes 
are accessible and can be controlled. This is an un- 
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realistic assumption in practice, as only a finite subset 
k = {fa, fa, • ■ •, fav} of the modes can be addressed si¬ 
multaneously. Moreover, as the estimation strategy is to 
be optimized at a fixed investment of energy, it appears 
unwise to initially populate modes in the complementary 
subset -> k, which cannot be measured. We therefore as¬ 
sume that the preparation of nontrivial initial states is 
also limited to k. Hence, the input states that we consider 
are of the form | ip) = \ip) k \0)^ k . The unitary U(9), 
on the other hand, acts on all modes such that the fi¬ 
nal state is | ip) = U(9)\ip), but only the reduced state 
Pk = Tr-,fe | ip) (ip | is accessible for the estimation of 9 
(see Fig. 1). Thus, Eq. (5) provides an upper bound, 
i.e., X(p k (9)) < I(\ip(9))) for the precision with which 9 
can be estimated. However, a precise expression for the 
losses incurred by tracing can be established as we will 
show now. Let us expand the reduced state Pk(9) in pow¬ 
ers of 9, 

Pk(9) = P? + 9^ +9 2 ~tf +0(9 3 ), ( 6 ) 


where the leading order is p k ] = Tr^j. (|'i/> (0) )(^ (0) |) and 
the corrections are given by 

p£= Tr.fcd^X^H-l^oX^WI), (7) 

p { k = Tr^ fc (| V> (0) }(^ (2, | + \ip ( - 2> )(ip ( - 0 ' , \ + 1 1 /; (1) )(^ (1) |) . (8) 

As the modes are noninteracting, the free evolution t/ (0) 
is a local operation that leaves the vacuum invariant and 
we may write |# 0) ) = U^\ip) k \0)^ k = \ip w ) k 1Ofa fe . 
It then easily follows that /5{ 0) = | ip^kki V’ (0) l is a pure 
state. In a similar way we may expand pk(9 + d9) as 


p k (9 + d9 ) = p k (9) 


, 0 dpk(0) d9 2 d 2 p k {9) 3 

dd ^^ + ^^9^ +0{dd ] 


= Pk + d9~p + d0 2 p™ + 0(9, dO 3 ). (9) 


Inserting (6) and (9) into Eq. (1) and taking into account 
the unitarity of the transformation, we arrive at the cen¬ 
tral result 


l(~Pk(9)) = ^ k (^\p^\r\ + 0(9), (10) 

for which a detailed proof is given in the Appendix. We 
pause briefly to highlight the simplicity of Eq. (10): One 
only requires a single matrix element in the second-order 
expansion of Pk(9) in order to compute the QFI, to lead¬ 
ing order in 9, for any input state | ip) k - 

To compare Eq. (10) to the case where no information 
is lost to tracing, we decompose p ( k as in Eq. (8). Ex¬ 
ploiting the unitarity of the transformation, one obtains 
(see the Appendix) 

l('Pk(9)) = X(|^» - A^k^um^m 2 + 0(9) 

i^O 

= 1(\$)) - AM) k ) + 0(9), (11) 

where I(\ ip)) is given as in Eq. (5) and the sum runs over 
all basis states of the Fock space of —ik that are orthogonal 


to the vacuum state |0fa fc . Note that Eq. (11) suggests 
that in order to compute the QFI to leading-order in 9 
one requires only leading order corrections to the initial 
state, whereas Eq. (10) above implied that second-order 
corrections are needed. However, as is shown in detail 
in the Appendix, the requirement that the transforma¬ 
tion is unitary implies that the second-order correction 
in Eq. (10) can be expressed in terms of the linear order 
corrections to the state. 


C. Parameter estimation from Bogoliubov 
transformations 

For the remainder of this paper, we will work with 
Eq. (11) to investigate how tracing losses can be mini¬ 
mized and optimal scaling can be achieved. For this task, 
we will specialize our discussion to Bogoliubov transfor¬ 
mations [30], that is, unitary transformations on the Fock 
space that can be viewed as linear transformations of the 
mode operators, i.e., 

a m l— t a m = 'y ] ( a mn a n ~ fimn a ri) > (1^) 

n 

where unitarity places restrictions on the complex coeffi¬ 
cients otmn and /3 mn . Another way to view such transfor¬ 
mations is as unitaries that are generated by quadratic 
combinations of the mode operators. In the perturbative 
regime that we consider here, the Bogoliubov coefficients 
are expanded as 

0 + 02 + °( 03 ) - ( 13 ) 

Pmn=^ n 0 + ^J 2 +O(9 3 ), ( 14 ) 

where a^ n = 5 mn G n = S mn exp(icp n ) encodes the 
phases (j) n that are accumulated due to the free time evo¬ 
lution t/ (0) . The coefficients a mn and /? mn are responsi¬ 
ble for the shifts of single excitations between modes m 
and n or the creation and annihilation of pairs of parti¬ 
cles in these modes, respectively. Any initial state of the 
Fock space may simply be written in terms of linear com¬ 
binations of products of creation operators acting on the 
vacuum. The operators transform according to Eq. (12), 
whereas the vacuum transforms as 

| 0 >-> 10 ) +0(9 2 ), ( 15 ) 

P,Q 

where G*/W J * = and /3^ must be a Hilbert- 

Schmidt operator J2 P q PpJ < oo to guarantee unitarity 
(see, e.g., Ref. [31, pp. 93 and 101]). We now determine 
the tracing losses A tI (\ip) k ) and the scaling of I(\ ip)) 
with the average energy of the initial state in the above 
scenario. 

The linear order operation t/ (1) can be (uniquely) sep¬ 
arated into (all) terms that leave the particle content of 
k or —ifc invariant and a term that correlates both sets 
via the creation or exchange of excitations, i.e., 

U (1 > = 0 + U ( k 0) 0 U^l + U { k % . ( 16 ) 
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From Eq. (11) it is obvious that U k ^®U^ k does not con¬ 
tribute to A tr (| V’)fe)- The other two terms generally pro¬ 
vide nonzero contributions. However, to linear order, the 
operation U <3) k populates |0)^ fc with pairs of particles, 
whereas U k \ k may either create one particle each in k 
and -iA; or shift one excitation from k to ~^k. As a result, 
for the same fixed state | i)_, fc , which contains at most 
two excitations, ^ k (i \ | 0 )^ k and k (4> m U k (i \ U£^ k \ ^) 

cannot both be nonzero and hence contribute separately 
to the tracing loss. The latter term can be avoided by 
choosing 11/’ } k such that superpositions of states that dif¬ 
fer only by one excitation are excluded. The former term, 
on the other hand, cannot be avoided as it is independent 
of the state | ip) k - It can, however, straightforwardly be 
computed using the vacuum state of Eq. (15), for which 
k(0U(i\U^ k \0) k \0)^ k = 0. For this case we find (see 
the Appendix) 

Atr(W fc )> A tr (|0) fc ) = 2 £ |/W>| 2 . (17) 

p,q(£k 

Thus, by choosing our initial state | ip) k such that su¬ 
perpositions of states that differ only by one excitation 
are excluded, the inequality in Eq. (17) becomes a strict 
equality and the QFI can be precisely determined to lin¬ 
ear order in a straightforward way using Eq. (11). We 
will restrict our discussion to such states from now on. 


D. Optimality 


We now show that the optimally achievable precision 
in this setup is Heisenberg scaling and how it can be re¬ 
alized even with nonentangled states. To determine the 
optimal states, the exact values of the Bogoliubov coeffi¬ 
cients a mn and (3 mn have to be known. This makes the 
procedure highly task specific. Nonetheless, it is possi¬ 
ble to identify families of states that scale optimally with 
the input energy under some minimal assumptions about 
the Bogoliubov coefficients. First, consider a single-mode 
Fock state with n excitations. Using Eqs. (12)-(15), this 
state transforms to 


\n k ) i ^ G k | n k ) 

- nG k a k l | nk) + \ \J(n + l)(n + 2) G k /3 kk 
^2(ynG k a ( k ^\n-l k )\l p ) - Vn+lG k /3^ 

p^k 


,(1) * |n 
.(!)„, 


n- 


p^k 

+ 2 l n fc)|lp)|l?) + Gpfi, 

p,q^tk p=£k 

p^q 


i(l)» 

'PP 


+ 0(9 2 ). 


A s l a fcpl = \ a pk\’ one immediately obtains the 


T(l nk)) 


= 2 n(n + l)!^! 2 + 4 n ^(|a$ | 2 + | 2 ) 

+ Z(|0» + 0(0). ^ (19) 


Notice that, here, the coefficient a(V, which leaves the 
occupation number of the mode k unchanged, does not 
contribute to the QFI, but it may do so for superpositions 
of different particle numbers. However, when f3 kk ^ 0, 
which corresponds to single-mode squeezing transforma¬ 
tions, terms that shift the population of the mode k 
by two excitations yield QFI that scales quadratically 
with n, i.e., Heisenberg scaling. On the other hand, clas¬ 
sical states, which in the present context are all coherent 
states, yield QFI that scales linearly with the average 
particle number [24]. 

As an example for the application of Eq. (19), con¬ 
sider a scenario in quantum optics [32] where only a single 
mode labeled k is subject to weak single-mode squeez¬ 
ing, while all other modes are left invariant. The cor¬ 
responding Bogoliubov transformation is represented by 
the coefficients a kk = cosh(0) = 1 + \9 2 + 0(9 3 ) and 
f3 kk = sinh(0) = 9 + 0(9 2 ). The optimal state for the 
estimation of 9 1 is hence already \n k ) and the QFI is 

I(| n k )) = 2 n(n + 1) + 0(0). 

When the diagonal linear coefficients /3 kk vanish, as in 
the transformation of field modes of nonuniformly accel¬ 
erating rigid cavities [25, 26], excitations in a single mode 
are not sufficient to obtain Heisenberg scaling and one 
requires initial states with at least two occupied modes. 
For instance, when the modes k and k' can be controlled 
and P kk = /3 k , k , = 0, the QFI for the state \n k )\m k ’) is 
(see the Appendix) 

l(\n k )\m k ')) = 8mn(\a^,\ 2 + \P kk ,\ 2 ) +4n^(|a^| 2 

P^ik 

+ 1/3^ I 2 ) +4m^ {\ a p k ' 1 2 + \Pp k '\ 2 ) + Z(|0» + 0(0). 

p** (20) 

The factor rrm allows the QFI to scale optimally with 
the average number of excitations for the estimation 
of beam-splitting (a^V,) and two-mode squeezing (/3 j f 1 . ) ,) 
terms. Indeed, inspection of the first-order expansion 
of |rife)|mfe/) reveals that Heisenberg scaling is the ulti¬ 
mate achievable precision scaling within our approach 
(see the Appendix). For the explicit example of a pure 
two-mode squeezing transformation, which may be re¬ 
alized in a variety of physical systems, e.g., in super¬ 
conducting circuitry [33], the Bogoliubov coefficients are 
given by a kk = a kV = cosh (0) = 1 + \9 2 + O(0 3 ) 
and /3 kk i = /3 k > k = sinh(0) = 9 + O(0 2 ). We hence re¬ 
cover the well known result that in this case the state 
| n k )\n k ') is optimal and the QFI from Eq. (20) yields 
l(\nk)\nk')) = 8 n(n + 1) + 0(9). 

Whereas the scaling with respect to the average en¬ 
ergy cannot be better than quadratic, even for arbi¬ 
trary Gaussian transformations beyond the examples 
presented so far, one can still improve the constant pref¬ 
actor in this scaling. For example, assuming fi k l = 0 
and a kk'i^kk' ^ O’ P ar bi c l e numbers n and n ± 2 
in the state (\n k )\n k ') + |n fc )|n-2 fe /) + |?r fc )|n+2 fe ^))/v / 3 
guarantee additional optimally scaling terms in the QFI 
while avoiding losses originating from |(-!/> (0) | ^ (1 >)| 2 . The 
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latter is nonzero whenever | ifi) contains superpositions 
of Fock states that can be converted into each other by 
changes of at most two excitations, such as \nk)\n+lk') 
and \ n+l k )\n k '), or |n fc )|n fc ') and | n+l k )\n+l k ')- 

Notice that in all of the examples considered so far the 
optimal states are not entangled, but can be regarded as 
highly squeezed and hence nonclassical. The final state 
may become entangled due to the transformation, but, as 
seen from (18) and (19), optimal scaling may arise from 
terms (here |ra—2*,)) that do not produce entangle¬ 
ment. Indeed, entanglement of the initial state is not 
necessary for optimality (a feature that was also noted in 
a related context in [34]), although one may find entan¬ 
gled states that do admit Heisenberg scaling in precision 
such as the state (| 7i+U)| n— 1*,') +1 n—1*,}] n+ lfc')) /a/ 2, 
which also has minimal tracing loss. For the entangled 
state above, homodyne measurements would be required 
for optimal scaling, whereas the measurements for the 
previous examples need to be number resolving. The 
choice of input state and measurement is thus also a mat¬ 
ter of what is practical in a given experimental setup. 

IV. CONCLUSION 

We have shown that the ultimate limits for estimating 
small parameters encoded in the linear transformations 
of the mode operators describing a set of noninteracting 
harmonic oscillators scale inversely proportional to the 
average input energy of the modes. Using a perturbative 


approach for these Bogoliubov transformations, applica¬ 
ble when the parameter of interest is small, e.g., as in 
the case of nonuniformly moving cavities, we have de¬ 
rived analytical formulas for the QFI, which apply to 
initial overall pure states of all modes. We have pro¬ 
vided a lower bound on the losses when some of these 
modes cannot be controlled. Finally, we have identified 
the class of states that yield optimally scaling precision 
while exhibiting minimal tracing losses. This provides a 
significant advancement beyond previous analysis of this 
problem [24], which was limited to pure Gaussian states. 
An investigation of mixed initial states, such as the mixed 
Gaussian states that were investigated recently [35] , goes 
beyond the scope of the current paper, but is certainly 
of interest. 

Our results open up the possibility to explore optimal¬ 
ity in a range of specific applications, where information 
about the Bogoliubov coefficients is available. Exam¬ 
ples include field modes in non-uniformly moving cavi¬ 
ties [25, 26], analog gravity phenomena [14], and effects 
in superconducting materials [15, 17, 36]. 
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Appendix: Explicit calculations 

In this appendix we present explicit proofs of some key 
results presented in the main text. In particular, we give 
a detailed derivation of the QFI for the estimation of 
small parameters 0 , encoded in Bogoliubov transforma¬ 
tions, and we discuss the optimality of Heisenberg scal¬ 
ing. 


A.l. Reduced state QFI 


We start with the QFI X(p k (d)) for a mixed state 
p k (8), which we will show to be given by 

l(p k (d)) = -A k (^\p^ \^°\ + 0(8) , (A.l) 

where p k (8) = Ti% fe | ^) (tjj |, |?/>) = U(8)\ip), and U(8) 
have small-parameter expansions of the form 

Pk(8) = pp + 8 + 8 2 p< 2) + G(8 3 ) , (A.2) 

11£) = |Vi (0) > + 0 P \V 2) ) + 0(8 3 ), (A.3) 

U(8) = U m + 8 U (1) + 8 2 U {2) + 0(8 3 ), (A.4) 

with |i/im) = \i/>) k ]0}^ = | '*/i (0> )fc 1The QFI 

of the state pk may be expressed via the Bures distance 
as 


Z(pk(0)) 


1 - 

lim 8- 

dQ—t 0 


^T(p k (0),Pk(0 + d8)) 

dd 2 


(A.5) 


where the Uhlmann fidelity is given by 

T(p k (8),p k (8 + d8)) (A.6) 

= (Tr ^JVh(8)p k (0 + d8)^p^)^ 2 . (A.7) 


With Eq. (A.2) at hand, the expansion in terms of d8 can 
be seen to be of the form 


p k (8 + dd) = p k (8) + dd 


dp k (0) 

88 


dd 2 d 2 p k (8) 
+ 2 88 2 


0(dd 3 ) 


= ~PT + de ~P { k + d^P^k + 0(8, d8 3 ). 

(A.8) 


Using Eqs. (A.2) and (A.8), we can then write 


\/ p k (8) Pk (8 + dd)\/ p k (8) 


I ~(0) ( ~(0) , 

= \]P k (Pk + 


de-tf + de 2 ~ P ^ 


+ 0(8, dd 3 ) 


= + d8 k (^\p£\r\ 

+ dd 2 k (^\pi 2) \^°\) + 0(8, d8 3 ), 


(A.9) 


where we have made use of the fact that p { k = 
I ' i / ;(0) )fcfe(V’ (0) l is pure. Now we can take a closer look at 
the term linear in dd, i.e., 

k(^\p?W°\ 

l(^ (1) l^)(fe(V’ (0) l ->fc(01) +H.c.)| i)^ fc |V’ (0) ) fe 

i 

= k (tp w \^ k (o\u^\ip) + (tp \ u (I)t |V’ (0) ) fe |o)_, fc 
= (V’l U (0)t t/ (1) IV) + (V>|U (1)t U (0) IV’) , (A.10) 


where we have used |'0 to) )fc 1 0 )_, fc = |i/i (0) ) = U m \'if>). 
On the other hand, the unitarity of the transformation 
requires 


1 = tfU = t/ (0)t t/ (0) + 8(U (0) m^ + t/ (1)t /7 (0) ) 

+ 6l 2 ([/ (0)t U (2, + £/ (2)t t/ (0) + U (1,t £/ (1) ) + 0 ( 8 3 ). 

^ (A.ll) 

Since t/ (0)t I/ (0) = 1, it follows that all other terms on the 
right-hand side of Eq. (A.ll) must vanish. We hence find 

k(r o) \p { k ) W°\= 0- (A.12) 

Reinserting this into Eq. (A.9), we obtain the Uhlmann 
fidelity 

T(p k (8),p k (8 + d8)) = 1 + dd 2 k (^\Pk } I 4> ( °\ 

+ 0(8, dd 3 ), (A.13) 

which in turn yields the QFI of Eq. (A.l). 


a 
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A.2. Relation to pure state QFI 


A.3. Tracing Loss 


Next we show that the QFI of the mixed state p k (0), 
arising from tracing out some of the modes from the pure 
state \ip), that is, p k = Tr^ fc \4>)('ip\, may be related to 
the QFI of the latter pure state via the expression 

z(pk( 0 )) = - iJ 2 \^ w U(i\u^m 2 + o( 0 ) 

i^O 

= l(\t)) - AM) k ) + 0(0), (A.14) 


where 

1(\4>)) = 4((# 1 >|$< 1 >) - |(Vi (0) |Vi (1) )| 2 ) +0(0). 

(A.15) 

Starting from Eqs. (A.l) and noting that p k ] can be writ¬ 
ten as 

pf = T>^(ivi (o >)<# 2, i + i# 2 ) x# o, i + , 

(A.16) 


we compute 

k(^ (0) \ PT\^ ( °\ = fc(^ (0) l Tr -fc |V’ (0) )(V’ (2) | 

+ |^ f2) )(Vi (0) | + |# 1) )(# 1) |] \^°\ 

= k( X (0) l * I [ (l 1 0 G (2,t + H. c.) 


U m 


U (1)t 


\Zk\^°\ 


= (V>|t/ (2)t f/ (0) + {7 (0)t G (2) \4>) 

+Y,\k(4> io) \^(i\ ua) m 2 ' 


(A.17) 


For the first term we can use the unitarity of the trans¬ 
formation from Eq. (A.ll) such that 

(V>|t/ (2)t G (0) + t/ (0)t t/ (2) 1 1/}) = - (ip\u (1 "u a) \ip) 

= -(4>w\ft 1 '>) . (A.18) 

In the remaining part, the vacuum state |0)^ fc may be 
extracted from the sum such that 

EI ^ (0, l ^(i\u (1) | V’)! 2 = I X# 0) l . fc <01 g u) \4>)\ 2 

i 

+ E \ k (i’ m \^(i\ uw \*l ’)\ 2 

ijto 

= \(r o) i x (i, )i 2 + E i *w (0) i • (a. 19 ) 

0 


We now specialize to the case where U(0) is a Bogoli- 
ubov transformation that is represented as a linear map 
of the mode operators, that is, 

a m l— t a m = ^ ' ( a mn a n ~ Pmn a ri) ) (A.20) 

n 

where the complex coefficients a mn and f3 mn can be ex¬ 
panded as 

= aZ + aZ 0 + a« d 2 + 0(0 *), (A.21) 

Pmn=PZ0 + ^n^ 2 +O(0 3 ). (A.22) 

In such a scenario, the vacuum state can be shown (see, 
e.g., Ref. [31, p. 100]) to transform as 

|0) —t |0) -0\Y. G lPm a l a \ l°> +°(° 2 )’ ( A - 23 ) 

p,q 

where G'/W 5 * = G*/?^)* due to unitarity (see, e.g., 
Ref. [31, pp. 93]). To evaluate the tracing loss for the 
vacuum, for which this loss is minimal, note that the cor¬ 
rection U w can be uniquely split up into terms 
and U^® U^ k , which act trivially on the subspaces k or 
-i k, respectively, and the term U k \ k i which correlates 
them, that is, 

U w = ® + U ( k 0) ® U™ + U ( Z k ■ (A.24) 

For the initial vacuum state 1 ip) k = |0) fc , only the term 
U k ) ® contributes to the tracing loss, which we iden¬ 

tify from the vacuum transformation as 

K 0) ® ^=4E G ^T a >t■ ( A - 25 ) 

p,q(£k 

With this, we may write the tracing loss for the initial 
vacuum state as 

Atr(|0) fe ) = 4El fc (0U fc <i|G-r®C^|0) fc |0K fc | 2 

= 2 E i§E g W^°i w 44 i°ui 2 » 

p',q'(fck p,q(£k 

(A.26) 

where we note that, as opposed to the sum over i, the sum 
over all p' and q' counts every basis vector twice and a 
factor of ^ therefore has to be included in the conversion. 
Since unitarity requires G*/3 pq = G*/3 qp (see Ref. [31, 
pp. 93]) and (Q\ a p m q 'a\,a\\0) = 5 pp '5 qq ' + 6 pq >5 qp > we 
arrive at 


Combining (A.17)-(A.19) and comparing with 
Eq. (A.15), we then immediately arrive at the re¬ 
sult of Eq. (A.14). □ 


Atr(|0) fe ) = 2 E l/^T ^ A M)k)- (A-27) 

p,q(£k 






A.4. Transformation & QFI for Example States 


A.5. Optimal scaling 


We now give the explicit expression for the state trans¬ 
formation of \n k )\m k >) and use the result to calculate 
the corresponding QFI. From Eqs. (A.20)-(A.23) we can 
write 

K)K'> -> G n k G%\n k )\m k ')+eGZG%' 

+ nG k a kk \ n k)\m k ') - \\/n(n-l)G* k \n-2 k )\m k *) 

- |v/(n+l)(n+2) G* k \n+2 k )\m k ') 


+ mG* k ,aft k , \n k )\m k i) - \\Jnn(m-l)G k ,p$ k , \n k )\m-2 k >) 

- \\J (m+l)(m+2) G* k , \n k )\m+2 k >) 

+ \Jm(n + 1) G k ,a$ k \n+l k )\m-l k >) 

+ \Jn(m + 1) Gla^l, \n-l k )\m+l k <) 

- sfrrvnG k /3 kk , \n-l k )\m-l k -) 

- \J (n+l)(m+l) \n+l k )\m+l k ') 

+ {y™ G k a kl \n-l k )\m k .)\l p ) 

p^k,k' 

+ VmG^all \n k )\m — l k '}\l p ) 

- Vn+IGIP^* |n+lfc)|mfc/)|l p ) 

- y/m+lG* k ,p™? |n fe )|m+l fe /)|l p )j 


l E 

p,q^k,k' 

Pt^Q 


G qPpq* W)\m k i) |lp)|lq) 


- I E G ;C* K>K'> |2p)] +0(e 2 ), (A.28) 

p^k,k' 

where + G n a^ = 0 and G* m ^ n = G* n ^ due 

to unitarity. For the initial state \ip) = \n k )\m k ’) we find 

(■0C 1 >|^< 1 >) = |n(n + 1) | 2 + \m{m + 1)\{)EI, 


+ n 2 \a kk \ 2 + m 2 \a ( i ) |2 


T'fe' 


2mn(\^ (1) |2 


n 2^ 0 a: 

p^k 

1£ K, ■ 


+ i a 


pfc I 


mE (l« 

p^k' 


2 

(i) |2 


pfc 7 1 


k k ' 

(1) I 2 ) 

(!) | 2 \ 


+ l/3pfc' 


)|2 


(A.29) 


and |($ (0) | V’ (1) )| 2 = ^l^fcfcl 2 + «i 2 |a 




(A.30) 

We hence obtain the QFI for the state \n k )\m k >) to be 


( 1 ) |2 


Z(\n k )\m k >)) = 2n{n+l)\/3^\ + 2m(m + l)\P£, k , 


- 8 mn (| c 


I 2 . 

■fcfe'l 


I/?, 




-4nE(l a 


+ 4m E d c 


p^k' 


; (1 >l 2 

pk 1 


p+k 

iC'i 2 ) +2 £i« 


pfc l A | i ) 

0(9). 
(A.31) 


( i )|2 
P9 1 


Finally, we show that, within the perturbative ap¬ 
proach we pursue here, no better scaling of the QFI is 
possible. To see this, note that Heisenberg scaling arises 
from the linear order correction terms in Eq. (A.28) that 
are proportional to square roots of at most quadratic 
combinations of m and n. Each of these terms cor¬ 
responds to either a shift of excitations between the 
modes k and k ', indicated by coefficients a5v,, or the cre¬ 
ation (annihilation) of a pair of particles in these modes, 
corresponding to the coefficients (3 kk , (3 kk >, and /3w?,. In 
other words, terms in | ^ (1) ) that scale as n are the re¬ 
sult of changing the occupation number of two highly 
excited modes by one excitation each or the occupation 
of one such mode twice. To linear order in 6, no more 
than these two excitations may be changed. No initial 
state may hence have a linear correction | ip w ) that grows 
faster than linearly in n. 


This still leaves the option to consider superpositions of 
states with different occupation numbers, such that the 
average number of excitations remains fixed, for instance, 
the state (\n k )\n k ') + \n k )\n-2 k >) + \n k )\n+2 k >))/\/3. 
The cross terms appearing in | may increase 
(or decrease) the overall QFI. However, to linear or¬ 
der in 6 the corrections to every fixed-excitation state 
in the superposition are nonorthogonal only to the lin¬ 
ear corrections to states that differ by no more than 
two excitations. For example, the linear corrections to 
both \n k )\n k ') and \n k )\n—2 k >) contain a term propor¬ 
tional to \n — lfc)|n — lfc'), whereas the linear correction 
of \n k )\n+2 k i) cannot contain such a term, as can eas¬ 
ily be verified from Eq. (A.28). The possible gain in the 
number of such cross terms that is obtained from adding 
more terms to the superposition is therefore at most lin¬ 
ear in the number of terms, which is compensated by the 
requirement of normalization. Consequently, the leading- 
order correction (in 9) to the QFI is at most quadratic 
in the average occupation number, as we have claimed. 


The higher-order corrections may feature higher pow¬ 
ers of n. However, the perturbative nature of the cal¬ 
culations means that the results presented here are only 
valid as long as the higher-order corrections are negligi¬ 
bly small with respect to the leading order. Therefore, 
the perturbative approach is not uniformly valid in the 
average occupation number. In a regime where n is large 
enough such that 0 2 I(|^>))/4 ss 1 the perturbative calcu¬ 
lation no longer yields a reliable result. We can therefore 
conclude that, within the perturbative regime, Heisen¬ 
berg scaling is indeed optimal. 













